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The effect of dimerization on the random antiferomagnetic Heisenberg chain with spin 1/2 
is studied by the density matrix renormahzation group method. The ground state energy, the 
energy gap distribution and the string order parameter are calculated. Using the finite size 
scaling analysis, the dimerization dependence of the these quantities are obtained. The ground 
state energy gain due to dimerization behaves as u"" with a > 2 where u denotes the degree of 
dimerization, suggesting the absence of spin-Peierls instability. It is explicitly shown that the 
string long range order survives even in the presence of randomness. The string order behaves 
as v?'^ with (3 ~ 0.37 in agreement with the recent prediction of real space renormalization 
group theory (/3 = (3 — \/5)/2 ~ 0.382). The physical picture of this behavior in this model is 
also discussed. 
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§1. Introduction 

In the recent studies of quantum many body problem, the ground state properties of the ran- 
dom quantum spin systems have been attracting a renewed interest owing to the occurence of 
exotic phases resulting from the interplay of quantum fluctuation and randomness. The real space 
renormalization group (RSRG) analysis of the 5=1/2 random antiferromagnetic Heisenberg chain 
has shown that the ground state of this model is the random singlet (RS) state in 
which the spins form singlets randomly with distant partners. This has been confirmed numerically 
by means of the density matrix renormalization group calculation.^ 

On the other hand, the effect of randomness on the spin gap state of one-dimensional quantum 
spin systems has been extensively studied theoretically and experimentallyi'@'i'i'0'lllllll'Ei'0'EHll'lll'© 
in the relation with the doping effect in spin ladders and spin-Peierls systems. Hyman et al%i have 
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applied the RSRG method to the spin-1/2 random dimerized antiferromagnetic Heisenberg chain 
and have shown that the dimerization is a relevant perturbation to the random singlet phase. They 
also argued that the ground state of this model is the random dimer (RD) phase in which the string 
long range order survives even in the presence of randomness. 

It is the purpose of the present work to investigate the ground state of the spin-1/2 random 
dimerized antiferromagnetic Heisenberg chain numerically using the density matrix renormalization 
group method. In the next section, we introduce the model Hamiltonian. The result of the 

numerical calculation is presented in §3. The last section is devoted to summary and discussion. 

§2. Model Hamiltonian 

The Hamiltionan of the spin 1 /2 random dimerized Heisenberg chain is given by, 

N 

n = Y,2JiSiS^+i, I 5i 1=1/2, (2.1) 

i=l 

with 

Ji = J{l + {-lf-^u + 5i), (2.2) 
where 6i takes random values. Here we assume the distribution 

P(,.) = ^ for -»V2< *.<,.,. 72. ^^^j 

otherwise. 

In the following, we define the energy unit by setting J = 1 and assume n > without loss of 
generality. The numerical calculation is done for W = 1. The bond which connects the i-th site 
and i + 1-th site is called the i-th bond. The bonds (sites) with even (odd) i are called even (odd) 
bonds (sites). 

The ground state of the undimerized counterpart of this model (n = 0) is the RS phase as verified 
by the RSRG method for the distribution function of the effective bond strengthMHB' In this phase, 
the effective energy scale and spatial scale L of a singlet pair is related as F = Infi oc —\fL. 
This implies that the spins form a totally singlet state in a segment of average length L at the 
energy scale ri. The size of each singlet segment grows as the energy scale is lowered and finally 
all spins form a totally singlet RS state. This calculation is extended to the dimerizied case by 
Hyman et In this case, the similar description holds down to the energy scale T ~ 1/u. At 
this energy scale, the average spatial scale of each singlet pair is Xjv? . As the decimation further 
proceeds, however, the odd bonds are renormalized to zero and the ground state is decribed as an 
assembly of RS segments whose average length is proportional to 1/ii^. This phase is called the 
RD phase. The spin-spin correlation length is finite and proportional to 1/n^ but the energy gap is 
distributed down to zero. Therefore this phase is a Griffith phase. Hyman et a& also argued that 
the string order defined for the regular bond-alternating chain§l) remains finite in the presence of 



2 



randomness. They claimed that the string order is proportional to v?^ with /3 = 2Ef for small u. 
Quite recently, however, Hyman and Yan^ have given a different estimation as /? = (3 — \/b)/2. 

§3. Numerical Results 

3.1 Average Energy Gap 

In Fig. 1, the average of the logarithm of the energy gap In A is plotted against \fN for n = 
and < 60 using approximately 250 samples. Here the bar denotes the sample average. For all 
these samples, the number m of the states kept in each DMRG step is 100 and only infinite size 
iterations are carried out. For the numerical caluculation, the samples are grouped into several 
subgroups. The error bars are estimated from the fluctuation among these subgruoups. For = 0, 
the average In A is well fitted by In A ~ 1.77 — 0.8\/iV as shown in Fig. 1. This also implies that 
the average In A is calculated accurately enough with m = 100. Considering that the gap tends to 
increase with u, the accuracy for the cases n > is expected to be even better than the undimerized 
case u = 0. Therefore we take m = 100 for the calculation of In A for finite u. Hyman et al% have 
shown that In A should behave as l/u for finite u and correlation length is proportional to l/v?. 
Therefore we may assume the following finite size scaling (FSS) formula. 

ln(A/Ao) = i/(nViV) with In Aq = 1.77. (3.1) 

The FSS plot is shown in Fig. 2. The data collapse onto a single curve fairly well confirming the 
conclusion of Hyman et a/.i-' We therefore analyse our numerical data assuming the FSS plot with 
the argument uN^/'^ in the following. 

3.2 Ground State Energy 

The average ground state energy per site Eg/N is also calculated for < 60 using approximately 
1200 samples for various values of u. The average magnetic energy gain per spin AEq/N due to 
dimerization is defined by {Eg{u) — Eg{u = 0))/N. To check the accuracy, we have carried out the 
calculation with m = 100 and 60 for 50 samples. The m-dependence was negligible compared to 
AEg even for u = 0.01. Therefore, we have taken m = 60 for the calculation of the ground state 
energy. Assuming the FSS formula, 

AE^/N = N-''/^g{u^), (3.2) 

we find good fit to a single curve for a ~ 2.46 it 0.10. Fig. 3(a) shows a FSS plot with a = 2.46. 
This indicates that AEg/N behaves as -u2.46±o.io small u in contrast to the regular dimerized 
chainii) in which AEg/N cx n^/^. Especially, it is clear that the FSS plot with a = 2 does not 
collapse onto a single curve as shown in Fig. 3(b) which ensures that a is definitely larger than 2. 
This implies that the spin-Peierls instahlity for infinitesimal u is unlikely in the random Heisenberg 
chain because the cost of the lattice deformaton energy is proportional to v?'. 
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3.3 String Order Parameter 



The string order parameter Og^p is defined for this system by extending the definition for the 
regular system as foHows, 



Ostr= lim 0^^,{l-N), (3.3) 

where Og^j.(/; A^) is the string correlation function in the chain of length A*" which is defined only 
for odd I as, 

ivr \ >N, (3.4) 

k=2i+l J 

where < ... >Ar denotes the ground state average for each sample of length A^. Although Hyman et 
al^ defined the string order parameter in somewhat different way, both definitions are equivalent 
in the thermodynamic limit. We have chosen this expression because the numerical error bars 
estimated as described in §3.1 are less pronounced. The sites k = 2i + 1 and 2i + I + 1 are fixed 
as distant as possible from the boundary for each / to reduce the boundary effect. The average is 
taken over 400 samples with A^ < 60 keeping 100 states in each step. Comparing the cases with 
m = 60 and 100, we have found only slight difference. Therefore we use the data with m = 100 in 
the following analysis. 

Figure 4 shows Os^j.(/;A^) as a function of / for relatively large values of u{= 0.1,0.2 and 0.3). 
It should be noted that 0^^^{N — 1; A^) = 1 because Y^k=i •^'1 = iii the ground state. For these 
values of u, it is clearly seen that the string long range correlation remains finite. To determine 



the precise u-dependence of Ogtr for smaller values of u, we resort to the FSS analysis. The string 



order parameter for the finite system Ostr(-^) is defined by Og^j-(A^/2 — 1; A^). We assume the FSS 



formula for Ogtr(-^) 

O^(iV) = N-l^hiu^). (3.5) 
For f3 = 0.37ib0.03, the data fall on a single curve within the error bars. In Fig. 5(a), the FSS plot is 



shown with /3 = 0.37. This indicates that Og^p is proportional to yO-''4±o.06 ^^le thermodymamic 
limit. This is in good agreement with the new estimation of /3 = (3 — \/5)/2 in ref. |^. Figure 5(b) 
shows the FSS plot with /3 = 2 as predicted by Hyman et a/.i^ at first. All curves do not collapse 
at all. Thus the possiblity /? = 2 is totally excluded within our numerical data. 



This exponent /3 is related to the exponent r/ of the size dependence of Ogtr(-^) n = in 
the following way. The RD ground state can be regarded as an assembly of almost decoupled RS 
segments. The left ends of each segment is the odd sites and right ends, the even sites. (It is assumed 
that the sites are arranged in increasing order of i from left to right.) The string correlation between 
the spins on both ends of a segment is perfect because J2i '^i — within each RS segment. Therefore 
this contribution is proportional to the square of the density of these boundary spins. Hyman et 
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al!& estimated (3 = 2 counting the contribution only from these bounadry spins. However, the finite 
string correlation can remain between the boundary spins and those inside the RS segments. This 
contribution to the string correlation between the 2i + f -th spin and 2j-th spin for large j — i >> 1 
is estimated as follows. 



Ostr(^;iV) = <exp<|ivr| ^ 5| + 

yk>2i+l;k&I k€l+l,..,J-l 



+ E Sl]\>^. (3.6) 

k<2j;k&J J ) 

where 2j = 2i + I + 1 and the capital indices distinguish the RS segments. It is assumed that the 
2i + 1-th and 2j-th spins belong to the different RS segments denoted by / and J, respectively. 
This expression is approximated as 



~ < exp 1 


1 E Sk \ >N 




[ fc>2i+l;feG/ J 


X < exp 1 


ivr E sA>r, 




[ fce/+i,. J 


X < exp 1 


ivr E sA>N. 




[ k<2j;keJ J 



(3.7) 

because the correlation between the segments is weak. The second average is estimated as 

< exp J ivr E ^k \ >N^ 1, (3.8) 

because these spins form an almost complete singlet state in each segment. If we assume that the 
string correlation between the boundary spin and the inner spins in a single RS segment decay by 
the power law as 

(3.9) 



< exp 1 


[i- E si] 

[ k>2i+l;kel J 


1 >7V 


< exp 1 


[i- E ^^1 

[ k<2j;keJ J 





^27-1^ (3-10) 

where d2i (c?2j-i) is the distance between the 2i-th (2j — 1-th) spin and the right (left) boundary 
spin of the /-th ( J-th) segment. Here the bars should be understood as the random average within 
each segment. Thus we find 



Ostr(Z;iV)~d27^27-i- (3-11) 
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The bar on the rhs represents the average over the whole sample. The averge size of the RS segment 
is proportional to Xjv} . Therefore the average string order Og^j- is proportional to u~'^'^ which yield 
/3 = 2??. 

In order to estimate the value of r/, we have measured the string correlation in the finite length 
undimerized random Heisenberg model. If the above picture holds, Og^j.(A^/2 — 1;A^) can be 
estimated between the + 1-th site and 3A^/4-th site as follows 



OstrW = OgtrW2-l,iV) 





r 3Af/4 "1 




< exp 1 


i^ E >^ 






[ fc=Af/4+l J 






r N N/i 




< exp 1 




E sA>M 




[ fc=l k=l 


fc=3Ar/4+l J 



= < exp hvr E >^< exp <^ ivr E >% 

y k=l J [ fc=3Af/4+l J 

~ (iV/4)-^(Af/4)-'^. (3.12) 

Here it is assumed that Y^k=i ^l. — ^ because the ground state is singlet. Thus we find, 

Ogtr(iV) « N'^^. (3.13) 

From the log-log plot in Fig. 6, we find 2ri ~ 0.4 it 0.05 which is close to the estimation /?(= 2rj) ~ 
0.37 determined from the u-dependence of Ogtr- 

§4. Summary and Discussion 

The ground state properties of the spin-1/2 random dimerized antiferromagnetic Heisenberg 
chain is studied using the DMRG method. It is verified that the average of the logarithm of the 
energy gap is inversely proportional to the degree of dimerization in agreement with the real space 
renormalization group theory. 

The energy gain per spin is found to be proportional to with a > 2. Therefore the spin-Peierls 
instability does not take place in the present model. This mechanism of the absence of spin- 
Peierls instability should be distinguished from the topological argument by Hyman et al% which 
excludes the possibility of spontanuously dimerized ground state in purely one-dimensional random 
spin system. In the spin-Peierls system, the phonons should be treated three dimensionally and the 
confinement of domain walls does not take place. It should be also emphasized that this argument 
concerns the spin-Peierls instability in the limit of weak spin-phonon coupling and does not exclude 
the possiblity of the spin-Peierls state for strong enough spin-phonon coupling. Actually, for weak 
disorder, it is reasonable to expect that energy gain AE/N behaves as n^/^ as far as u >> W. 
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Therefore, the spin-Peierls state would be stable against disorder, if the lattice distortion mq in the 
absence of randomness satisfies this condition uq >> W. 

The string order behaves as u^^ with (3 ~ 0.37 in agreement with the prediction of Hyman and 
YangS It is suggested that the discrepancy from the earlier estimation by Hyman et is due 
to the string correlation between the inner spins in different RS segment which was not taken into 
account in ref. ^. This picture is consistent with the behavior of the string correlation in the 
undimerized random Heisenberg chain. 

In this work, we have fixed the strength of randomness W = 1. Nevertheless, we expect our 
conclusion is universal at least qualitatively as long as W » u. It should be noted, however, the 
physics can be totally different for other kind of randomness such as the case of topological disorder 



studied in refs. 15 and 16 



The numerical calculations have been performed using the FACOM VPP500 at the Supercom- 
puter Center, Institute for Solid State Physics, University of Tokyo. This work is supported by 
Grant-in- Aid for Scientific Research from the Ministry of Education, Science and Culture. 
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Fig. 1. The system size dependence of In A for m = 0. 



Fig. 2. The FSS plot of In A for u = 0.01 (o), 0.02 (•) and 0.03 (□). 



Fig. 3. The FSS plot of AEg/N for u = 0.01 (o), 0.02 (•) and 0.03 (□) with (a) a = 2.46 and (b) a = 2. 



Fig. 4. The spatial dependence of the string order OgtrC^; for u = 0.1 (o), 0.2 (•) and 0.3 (□). 



Fig. 5. The FSS plot of OgtrC-^) for w = 0-01 (o), 0.015 (□), 0.02 {o) and 0.03 (A) for (a) (3 = 0.37 and (b) 13 = 2. 



Fig. 6. The log- log plot of the string correlation OgtrC^) in the undimerized random Heisenberg chain(M = 0) 
plotted against N. The solid line is the fit by the power law OgtrC^) ~ N~'^'^ with 77 = 0.2. 
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